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We introduce a general class of generating functionals for the calculation 
of quantum-mechanical expectation values of arbitrary functionals of fluctu- 

in 

r^ . ating paths with fixed end points in configuration or momentum space. The 

o' 

y^ I generating functionals are calculated explicitly for harmonic oscillators with 

o ' time-dependent frequency, and used to derive a smearing formulas for correla- 

tion functions of polynomial and nonpolynomials functions of time-dependent 
positions and momenta. These formulas summarize the effect of thermal and 



^ ' quantum fluctuations, and serve to derive generalized Wick rules and Feyn- 



man diagrams for perturbation expansions of nonpolynomial interactions. 

I. INTRODUCTION 

A useful technique for describing compactly the properties of a quantum mechanical 
system is to define a suitable generating functional of some external source or current j{t). 
The desired properties are obtained from functional derivatives with respect to j{t). For 
example, the correlation functions and the quantum mechanical density matrix in one space 

1 



dimension x is determined by a generating functional which is a path integral in configuration 

space over all paths x{t) with fixed end points x(ta) = Xa,x(tb) = x^ |1], Chap. 2]: 

XbM 
{Xbh\xata)[3{t)] = I Px(t)exp|^^[x(t);j(t)]} , (1-1) 

where the exponent contains the classical action ^[a;(t)] plus a source term linear in x{t): 

tb 
A[x(t)-3(t)]= A[x(t)]+ j dtx(t)3{t) . (1.2) 

tb 



In this note we set up a useful alternative expression for the generating functional (|1 . 1| ) and 
a related one in momentum space. This alternative expression is obtained by extending the 
current jit) by singular sources proportional to 5{th — t) and 5{t — t^), and by reducing the 
path integral ( |1.1| ) with fixed end points in configuration space to one with vanishing end 
points. This will permit us to simplify considerably the calculation of quantum mechanical 
correlation functions. To see this simplification explicitly, consider a harmonic oscillator 
whose action reads 



A[x{t)] = dt 



ta 



ya;2(t) - y^V(t) 



:i.3) 



for which the generating functional can be calculated [0, Eq. (3.89)] as follows: 



, I UV+^l - / ^^ j ^^^ [(^b + ^g) COS a;(4 - tg) - 2XaXb] 

J i f ^^ Xasmu{tb -t) + Xbsmujjt - tg) 
X exp \t I '^^ '■ 77 r^; 3 v^i 

n J smujith — t„) 



ta 
tb t 



dt I dt-^^""(^^-^l^^""^f-^°^■(t)J^(o \ . (1.4) 



hMuj J J sinci;(tf, — t) 

The nonzero end points Xg and Xb make this expression quite involved. For vanishing end 
points, however, it simplifies to 



xr / XT / MuT 

{Oh\OQ[j{t)]-- 



I 2TTih sinuj (tb — tg) 

The observation which motivates the present paper rehes on replacing in the simple expres- 
sion ( |1.5|) the current j{t) by 

fit) = j{t) + MxJit - to) - Mxb6{h - t)], (1.6) 

where the delta functions are understood as 6{t — t^ + e) and 6{ti, — e — t) in the limit 
e — >^ 0. By performing some partial integrations, this replacement reproduces all terms in 
the complicated generating functional ( |1.4[ ), except for a rather trivial additional singular 
phase factor. The important relation is 

(Xfe h \Xata) b (t)] = (Xfe = tb I X„ = t„) [j (t) + MXa6{t - t J - MXb6{h - t)] 

xew[^{xl + xl)6{0)]. (1.7) 

In Section |I| we prove that the relation ( |1.7D holds for an arbitrary quantum-mechanical 
system whose Hamiltonian has the standard form 

Hoip,x,t) = ^ + Vix,t). (1.8) 

In Section |I| we calculate explicit amplitudes for a harmonic oscillator with arbitrary time- 
dependent frequency, and as an important application we derive in Section ^V\ from the 
new form of the generating functional a smearing formula for calculating expectation values 
of polynomial and nonpolynomial potentials functions. In particular, this result allows to 
calculate expectation values appearing in perturbation expansions for nonlinear interactions, 
as for example for the nonlinear a-model. In Section V we show that our smearing formula 



generalizes Wick rules and Feynman diagrams for harmonic expectation values from products 
of variables to mixtures of nonpolynomial functions and polynomials. In Section [V^, we 
finally specialze our generating functional simplifies to periodic paths. 

II. GENERATING FUNCTIONALS 

We begin by setting up phase space path integrals for generating functionals with fixed 
end points in either configuration or momentum space. The action contains additional 
currents k{t) and j{t) coupled linearily to momentum p{t) and position x(t). By extending 
the currents with singular (5-functions as in Eq. ( p..6|) , we reduce the path integrals with 
fixed end points to those with vanishing end points. Our procedure applies to arbitrary 
Hamiltonians HQ{p,x,t), with certain simplification resulting from a standard Hamiltonian 

A. General Phase Space Formulation 

Consider a quantum-mechanical particle coupled to a momentum and position source 
k{t) and j{t) with the classical Hamiltonian 

H{p, X, t) = Ho{p, X, t) - pkit) - xj{t), (2.1) 

where the corresponding action reads 

A[p{t),x{ty, k{t),j{t)] = J dt{p{t)xit) - H{p{t),x{t),t)}. (2.2) 

ta 

The total time evolution amplitude between fixed space points Xa and Xb is given by the 
path integral 



ix,h\xata)[kit),jit)] = J 27S exp|^^[p(t),x(t);fc(t),j(^)]} • (2-3) 

A Fourier transformation with respect to Xa and Xh produces the time evolution amphtude 

in momentum space 

+00 +00 

{Pbh\pata)[k{t),m] = J dXa J rfXfc C-^^^"^""^'"^")/'^ ( Xfc 4 | X„ t, ) [^ W, J W] • (2-4) 

—00 —00 

Here the initial and final momenta pa and pb are held fixed, so that the right-hand side may 
be written as the path integral 

iPbh\Pata)[k{t),jit)]= I ^^(^i|^exp|l^[p(t),x(t);Mi),jW]}- (2-5) 



Pa,ta 



We remark that both path integrals (|2.3|) and ( p. 51) are properly defined as continuum limits 
of ordinary integrals after a time-slicing procedure. Since end points of paths are fixed 
in coordinate and momentum space, respectively, the discretized expressions for the path 
integrals turn out to be slightly asymmetric in p(t) and x{t) [|l]. Chap. 2]. 

The time evolution amplitudes (|2.3|) and (|2.5|) with fixed end points can now be reduced 
to corresponding ones with vanishing end points. For this, we shift the current k{t) in ( p.l|) 
by a source term XbS(tb ~t)— XaS(t — ta) and observe that this produces by ( |2.2|) and ( p.5|) 
an overall phase factor: 

iPbh\Pata)[k{t) +Xb5{tb-t) -Xa5{t-ta),j{t)] 



TiVbXh - VaXa) 



{Pbtb\Pata)[k{t),3{t)]. (2.6) 



= exp 

By inverting the Fourier transformation ( p.4[ ), the configuration space amplitude ( p.3| ) is 
seen to satisfy 

{Xbtb\Xata)[k{t) +XbS(tb-t) -X'a6{t-ta),j{t)] 

= {Xb + x'b tb\Xa + x'a t, ) [fc(t), j(t)] , (2.7) 



where again the delta functions are understood as S(tb — e — t) and S(t — ta + e) in the hmit 
e — > 0. Because of this relation, the amplitude ( |2.3| ) can be reduced to a path integral with 
vanishing end points but additional (5-terms in the current k{t): 

(XbU I Xata)[kit),j{t)] = {xb = Oh\Xa = Ota)[k{t) + Xb5{U - t) - Xa5{t - ta),j{t)]. (2.8) 

A similar expression exists, if momentum end points are fixed in momentum space by adding 

Va^it — to) — Pb^itb — t) to the current j{t): 

{Pbtb\PataMt),j{t)] = {p, = Otb\Pa = OQ[k{t),j{t)+pJ{t - Q - Pb6{h " t)]. (2.9) 

We now explore the consequences of these two relations for the calculation of correlation 
functions. 

B. Correlation Functions 



The functional dependence of the time evolution amplitudes (|2.3|) and (|2.5|) on the 



currents k{t) and j{t) allows us to calculate expectation values of arbitrary functional 
F[p{t),x{t)] from the path integral 

{F[p{t),x{t)])[k{t),j{t)]l''\' = ^ \,,^^ , ,^ 

X / ^^^|^F[p(t),x(t)]exp|l^[p(t),x(t);A:(t),j(^)]} , (2-10) 

Va,ta 

where the variable v may be p or x. The usual correlation functions 

Vb,tb _ 



{P{tl) ■ ■ -piQxih) ■ ■■x{Un)mt),j{t)]Z'l 



{Vbtb\Vata)[k{t),j{t)] 

X J 2^^^ P^ti) ■ ■■P{tn)x{ti) ■ ■ ■x(t^) exp |^^[p(t),x(t); k{t),j{t)] \ , (2.11) 



are special cases of (|2.1CI| ), so we shall call the general expectation values ( |2.1(]| ) correlation 



junctionals. The sources k{t) and j(t) permit us to express ( |2.1(]| ) in terms of functional 
derivatives: 

'Ti 5 n 5 



Vb,tb 



{F[p{t\xm[Kt\m]z 



i6k{ty i5](t) 



{Vbh\Vata)[k{t),j{t)] 

. (2.12) 



{Vbtb\Vata)[k{t),j{t)] 

Recalling ( |2.8| ) and ( p.9|) , we shall rewrite the functionals {vbtb \ Vata)[k(t),j{t)] in a unified 
common way as follows 

{Vbtb\Vata)[k{t),j{t)] 

^^^^^^''^^h{v{ta)-Va)6{v{tb)-Vb)exp{^A[p{t),x{ty,k{t),j{t)]\ , (2.13) 



Wa=0,ta 

where the paths v{t) stand either for p{t) or for x{t). In each of these cases, the paths 
w{t) denote the conjugate variables x{t) or p{t), respectively. In this form, the path integral 
possesses the advantage that usual correlation functions (|2.11|) can be determined by path 
averages, in which intermediate and end points are treated on equal footing. Indeed, inserting 
delta functions according to 



{Pitl) ■ ■ ■P(t„)x(ti) ■ ■■xitm))[kit),j{t)]Z'l 

+ 00 +00 +00 +00 

- / dpi--- / dpn / dXi--- / dXmPl---PnXl---Xm 



{Vbtb\Vata)[k{t),j{t)] 



— oo — oo — oo 

X Y^^^§^Sip{t,)-p,)...6{p{t^)-p^) 

Va,ta 

x5(a;(ti)-si)---,5(x(tj-x„)exp|^^[p(t),x(t);/c(t),j(t)]^ , (2.14) 



we obtain with a similar reasoning 

{pitl) ■ ■ •p(t„)x(ti) ■ • ■x(t™))[A;(t), j(t)]:;::t 



+ 00 +00 +00 

1 



/ dpi--- / dpn / dXi--- / dXmPi---pnXi---Xr. 



{Vbtb\Vata)[k{t),j{t)] 

— OO — CXD — OO — CX) 

«)(,=0,t(, 

x5(x(ti) - xi) ■ ■ ■ 6{x{tJ - Xm)S{v{h) - t;fe) exp |^^[p(t), x(t); fc(t), j(i)] ^ • (2-15) 



C. Standard Hamiltonian 

The above formalism can be made more specific for the standard Hamiltonian ( |1.8D . 
Then the path integrals over the momentum paths p{t) in (|2.3|) and ( |2.5|) becomes harmonic 
and can be explicitly evaluated. The phase space integral (|2.3| ), for instance, reduces to the 
configuration space path integral 

{Xbtb\xata)[k{t),jit)]= I Vxit)expi^^A[x{t);k{t),j{t)]^ , (2.16) 

where the current k{t) couples linearily to the path momentum Mx{t) in the action 

A[x{ty, k{t),j{t)] = jdt [^x\t) - V{x{t),t) + x{t)j{t) + Mx{t)k{t) + yA^^W} • (2.17) 



ta 



A subsequent partial integration transforms the current k{t) to an effective coordinate cur- 
rent with an extra phase factor: 



{Xbtb\Xata)[k{t),j{t)] 

tb 



[Xbtb I XataW^jit) - Mkit)] exp { ^-j- 



Xbkh -Xaka + - / dtk'^{t) 



(2.18) 



In the next section we determine the generating functional {xb tb \ Xa ta) [0, j{t)] for a harmonic 
oscillator with arbitrary time-dependent frequency Q(t) and use ( |2.18| ) to construct the full 
generating functional {xbtb \ Xata)[k{t),j(t)]. 

8 



III. TIME-DEPENDENT HARMONIC OSCILLATOR 

Consider a standard Hamiltonian (|1.8|) with a harmonic potential containing an arbitrary 
time-dependent frequency: 

V{x,t) = ^n\t)x\ (3.1) 



The generating functionals ( |2.3| ) and ( |2.5| ) are then expressable in terms of two fundamen- 



tal solutions Da{t),Dh{t) of the corresponding classical equation of motion with particular 
boundary conditions 

K{t) Dait) = ; Daita) = , bait:) = 1 , (3.2) 

k{t)Db{t) = Q- Dbih) = , bbitb) = -I , (3.3) 

where K{t) denotes the operator 

Kit) = -d^-n\t). (3.4) 

Since the time derivative of the Wronski determinant 

W{t) = Da{t)bt{t) - ba{t)Db{t) (3.5) 

vanishes, we observe the identity 

Daih) = D,{ta) . (3.6) 

Note that a similar identity does not hold for the time derivatives of the two fundamental 
solutions Dait) and Db(t). Indeed, partially integrating the differential equation for Da(t) 
and taking into account (|3.2|)-(pl3D, we deduce 

9 



DbiQ + Mh) = -2 / dtnitw) Dait)Dhit) . 



(3.7) 



Let us now determine the time evolution amplitude ( p.l6|) in configuration space for a 
vanishing current k{t) as defined in (|2.16| ). We decompose the paths x(t) into the classical 
path xliit) and the quantum fluctuations 6x(t) around it: 



x{t) = xii{t) + 5x{t) . 
The classical path xli{t) solves the boundary value problem 

m 



(3.8) 



K{t)xUt) 



M 



■^clv «ij — ^a ; •^clv^f'/ — "^b ) 



(3.9) 



and the fluctuations Sx{t) vanish at the endpoints: 



6x{Q = 6x{U) = . 



(3.10) 



Inserting the decomposition ( p. 8] ) into the action ( |2.17| ), we observe that due to ( |3.8| ) and 
(p.9|) the total action decomposes into a classical part 



M r 



•^[a^cl(^); 0, j(^)] = TT XbXi^ih) - XaX^yita 



1 *' 

+ -ldtxi,{t)j{t), 

ta 



(3.11) 



and a fluctuation part, which is simply the classical action evaluated for the fluctuations 
6x{t) at j = 0: 



A[xity, 0, Jit)] = A[xUty, 0, Jit)] + A[6xit);0, 0] . 



(3.12) 



Inserting this into the original path integral ( p.l6| ), it factorizes into the product of a classical 
amplitude with the classical action ( |3.11| ) , and an additional fluctuation factor which is equal 
to the amplitude at vanishing end points: 



ixbtb\Xata)[0;jit)]=exp\-A[xi,it);0,jit)]\ixb = Otb\xa = Ota)[0,0]. (3.13) 



10 



A. Classical Action 

The classical action in the presence of currents can be expressed in terms of the solutions 
Dait), Db(t) of the time-dependent harmonic boundary value problems ( p.2[ ) and ( |3.3| ). First 
we decompose the solution of the boundary value problem ( |3.9| ) in the presence of into a 
homogeneous and an inhomogeneous contribution: 

xiM = ^ci{t) + Axi,{t) . (3.14) 

The homogeneous solution reads 

_ Dbit)Xa + Dait)Xb . 

while the inhomogeneous one is given by 

AxM = --ldt'G^^{t,t')j{t'), (3.16) 

ta 

where G^j{t,t') denotes the Green function of the classical equation of motion 

K{t)G';^{t,t') = 6{t-t'), (3.17) 

with Dirichlet boundary conditions 

6-^,(4,0 = G'J,(t,,t')=0. (3.18) 



/^From ( p.l7|) we deduce that the Green function G^j(t, t') solves the homogeneous differential 



L^ 



equation for t ^ t': 



K{t)G]^{t,t') = 0, (3.19) 



and that its first derivative dtG^j(t,t') is discontinuous at t = t': 

11 



lim 

40 



dtGUt,t')^ ^,, -dtGUt,t' 



t=t'+e 



n^ 



t=t'-e 



-1. 



The Green function itself is continuous around t = t': 



(3.20) 



lim 

40 



'^ii(^'^') ._„^. ~ '^ii(^'^' 



t=t'+e 



t=t'-e 



The solution of ( |3.18| )-( ]3.!21| ) is given by Wronski's famous expression 



^. (^^^,) ^ e(t - t')D,it)Dait') + ejt' - t)Da{t)D,{t') ^ ^,. ^^,^^^ 



Da{h) 



(3.21) 



(3.22) 



where 0(t — t') denotes the Heaviside function which vanishes for t <t' and is equal to unity 
for t > t'. Inserting (|3.14|) and ( p.l6|) we obtain for the classical action (|3.11|) 



A[T>,,{ty,o,j{t)] 



M 

2Da{h) 

tt 
dtxc\{t)j{t) 



Da{h)x^ - Db{ta)Xa - 2XaXb 

tb ib 

dt fdt'GUt,t')jit)jit') 



1 



2M 



where Xdit) and G^At,t') are given by (|3.15|) and ( p.22| ), respectively. 



(3.23) 



B. Fluctuation Factor 



Now we calculate the fluctuation factor in (|3.13|) . Recalling the path representation 
( p.l6| ) with the action (|2.17|) , we have to evaluate 



5x(,=0,ti, 

{xb = 0h\xa = 0ta)[0,0]= J V5xit)exp 

SXa=0,ta 



-- / dt5x(t)K(t)5x(t) 
2n J 

ta 



(3.24) 



To this end we decompose the fluctuations 6x{t) in ( ^.24] ) into eigenfunctions x„(t) of the 
operator K{t) of (|3.4|) with Dirichlet boundary conditions 



k{t) Xn{t) = Xn Xn{t) ] Xn{ta) = X„(tb) = 



(3.25) 



12 



which satisfy the orthonormahty and completeness relations 

tb 



/ dtXn{t) Xn'{t) = 6n,n' , (3.26) 

ta 

Y,Xn{t)Xn{t')=6{t-t'), (3.27) 



as follows: 



5xit) = J2cr.Xnit). (3.28) 



The path integral over all possible fluctuations 6x{t) in ( p.24|) amounts to a product of 
integrals over all expansion coefficients c„: 

5x{tt)=0 / +00 ^ 

f V5x{t) = J Ijl I dcn\ . (3.29) 

5x{ta) = I " -OO J 

The Jacobi determinant J of the transformation ( p.28| ) is an irrelevant constant. Applying 
( p.25| )-( p^ ), the path integral ( |3.24| ) is finally determined by 



(x, = Otfc|x, = OtJ[0,0]= ^ , (3.30) 

VDetir(t) 

where the determinant of the operator K{t) is equal to the product of its eigenvalues 

Det i^(t) = n^n. (3.31) 



C. Operator Determinant 



In order to calculate the operator determinant ( p. 311) it is advantageous to introduce a 
one-parameter family of operators 0,^ 

k^(t) = -d^,-g^\t), (3.32) 

13 



depending linearily on a coupling strength parameter g G [0,1], and coinciding with the 
original operator K{t) in ( |3.4| ) for g = I. It is possible to relate the operator determinant 
DetK^{t) to the fundamental solutions D^{t), D^{t), and to the Green function Gff(t,t') 
emerging from (|3.2|), { \i.'S\) and ( |3.17| ), ( |3.18| ). For this we substitute the operator K{t) 
by K^{t), and differentiate the (^-dependent version of the eigenvalue problem ( ^.251 ) with 
respect to g: 



k^{t) ^4^ - ^\t) xi{t) = ^ xi{t) + Xi ^"^'^^^ 



dg 



dg 



dg 



(3.33) 



Multiplying (|3.33|) with x^(t)/A^ and performing a summation over n plus an integration 
with respect to t, we obtain with ( |3.25|) , (|3.26| ) and (|3.31|) : 

rl i 

' -^ ^^^"'■' ' ^■^'"■' ^''"'- ^ (3.34) 



In Det K^t) = - dt Q\t) Gff{t, t) . 
In the last step we have used the spectral decomposition of the Green function 



^j'j [^^^ ) - 2^ Tg 



An 



(3.35) 



To solve the differential equation ( p.34|) , we differentiate the boundary value equation ( p^ 
for D^(t) with respect to g, and obtain the inhomogeneous initial value problem 



k^{t)^^ = n'it)Di{t); 



dDi{t) 



d dDi{t) 



t=ta 



dt dg 



0. 



(3.36) 



t=ta 



dg w av., Qg 

Generalizing ( p.22| ) from g = 1 to arbitrary values g G [0, 1], the solution of (|3.36|) is given 

by 



-lnD^,{U) = ~JdtQ\t)Gff{t,t). 



(3.37) 



This shows that (|3.34]) is solved by 



14 



DetK^{t) = CD^^{U) 



(3.38) 



where C denotes some constant. Due to this result, the ratio of two fluctuation factors with 
two different parameters gi and (72 can be rewritten as 



{Xb = 0h\Xa = 0ta)[0,0] 



91 



{Xb = 0h\Xa = 0ta)[0,0]3^ \ Di' (h) ' 



Di\U 



(3.39) 



This serves to determine the fluctuation factor of the initial time-dependent harmonic os- 
cillator at (71 = 1 in terms of the fluctuation factor of the free particle g2 = 0. The latter is 
well-known and may be calculated explicitly, for instance, via time-slicing |]1], Chap. 2] as 



Xb = Otb\Xa = Ota)[0,0] 



92=0 



M 



2Tcih{th - to) 



(3.40) 



Since the obvious solution of ( p.2|) at (72 = reads D^^ °(tf,) = t^ — ta, we obtain the famous 
Gelfand-Yaglom formula for Dirichlet boundary conditions H| : 



Xb = Otb\Xa = Ota)[0,0] 



M 



2TTihDa{tb 



(3.41) 



Note that similar results can also be derived for periodic and antiperiodic boundary condi- 
tions m. 



D. Full Generating Functional 



Having obtained the generating functional {xbtb \ Xata)[0,j{t)] of the harmonic oscillator 
with arbitrary frequency with vanishing current k{t), we now make use of the relation (|2.18 



to derive the full generating functional {xbtb\xata)[k(t),j{t)]. The terms containing the 
current velocity k{t) can be turned into functionals of k{t) itself with the help of several 
partial integrations. These turn out to remove the extra phase factor in (|2.18| ). As a result. 



15 



the time evolution amplitude in the configuration representation is determined by a Van 
Vleck-Pauli-Morette type of formula |I|, Chap. 4] 



{Xbh \Xata)[k(t),j{t)] 






(3.42) 



with the action 



A{Xb,tb;Xa,ta)[k{t),j{t)] 



M 



Da{tb)xl - bb{ta)xl - 2XaXb 



tb 



1 



2Da{tb) 

ib tb 



+ Jdt [x,i{t)j{t)+p,i{t)k{t)] -^J dtj dt' 



j^G^it^'WW) 



+G%{t,t')j{t)k{t') + G^t,t')k{t)j{t') + MGUt,t')kit)k{t') 



(3.43) 



The homogeneous classical solution Xdit) is given in (|3.15|) , and Pdit) denotes the classical 
momentum Pci(^) = ^a^ci(^)- The Green function G^At,t') is given by (|3.22|) , while the 
others are 



^. . , ,n _ e(t - t')Dbit)Dait') + ejt' - t)Da{t)Db{t') _ 

Gjkit, t ) ^-^^ G,^{t , t) , 

r^ (, ,,^ _ e(t - t')Db{t)ba{t') + 9(f - t)b,{t)bb{t') _ 



(3.44) 
(3.45) 



By differentiating (|3.43|) functionally with respect to j and k, we see that the Green functions 
correspond to the correlation functions 



(X6|x(t)x(t0ka) = — GJ^.(t,t'), 






( Xb I xit) Pit') \xa)= ifiG^it, t') = ihGl^{t\ t), 
{xb\p{t)p{t')\xa) = inMGlj,{t,t') 



(3.46) 
(3.47) 
(3.48) 



with x{t) = x{t) — Xci(t) and p{t) = p{t) — Pc\{t). These results can be summarized by the 
mnemonic rule that the Green functions involving a momentum current k{t) once or twice 



16 



follow from G^jAt^t') by one or two time derivatives if the time derivatives of the Heaviside 
functions are neglected: 






dt 



dt' 



dt.dt' 



(3.49) 



A complete analogous expression to ( p.43| ) is found for the time evolution amplitude in the 
momentum representation. The Fourier transformation ( |2.4]) of (|3.42|) yields a Van Vleck- 
Pauli-Morette type of formula 



{pbtb\Pata)[k{t),j{t)] 



'Sttz/i ^— ^ exp^-A{pb,tb;Pa,ta)[k{t),j{t)] 



(3.50) 



where the action is the Legendre transform of (p.43| ) 



A{pb,tb;Pa,ta)[k{t),j(t)] = A{Xb,tb;Xa,ta)[k(t),j{t)] - PbXb + PaXa, 



(3.51) 



calculated for the conjugate variables 



_ dA{Xb,tb;Xa,ta)[k{t),j{t)] _ dA{Xb, tb;Xa,ta)[k{t),j{t)] 

Pb — <-, 1 Pa — 



dxh 



dXn 



(3.52) 



This brings (|3.51| ) to the form 



A{pb,tb;Pa,ta)[kit),j{t)] 



Da{tb) Da{tb)pl - Db{ta)pl " 2paPb 



2M[l + Da{tb)Db{ta)] 

tb tb 



+ j dt [xci{t)j{t)+p,i{t)k{t)] -l^j dt j dt' 

ta (■a Cq, 



-G^,(t,t')j(t)j(0 



+G%{t,t')j{t)k{t') + GUt,t')k{t)j{t') + MGUt,t')k{t)kit') 



where the classical solution now reads 



(3.53) 



Xcl{t) 



PalDajt) + Dbit)Da{tb)] + Pb[Da{t)Db{ta) - Db{t)] 
M[l + ba{tb)Db{ta)] 



(3.54) 



17 



and Pci(t) denotes the associated classical momentum pdit) = Mxci(t). The Green functions 
in ( |3.53| ) turn out to be 



GUt,t') = eit-t') 



,. [D,{t)Daih) + Da(t)][Da(t')b,{Q - D,{t')] 



+e(t' - 1) 



Da{h)[l+Da{h)D,{ta)] 
[Da{t)b,{ta) - D,{t)][D,{t')ba{t,) + Da{t')] 
Da{U)[l + ba{h)bb{ta)] 



Gjj [t ,t) , 



GUt,t') = Q{t-t') 



,, [D,{t)Da{U) + DamDa{t')D,{ta) " D,{t')] 



+e{t'-t) 



Daih)[l+Daih)D,iQ] 
[Da{t)bb{ta) - DbimMtlDajtb) + Dajt')] 
Da{U)[l + ba{h)b,{ta)] 



GUt,t') = Qit-t') 



+Q{t'-t) 



,. [b,{t)bM + b,mDa{t')bt{ta) - b,{t')] 



Daih)[l+Da{h)D,{ta)] 
[ba{t)bb{ta) - b^{t)]{b,{t')ba{U) + bait')] 

Da{h)[i + ba{h)b,{ta)] 



The relation to the correlation functions is similar to ( |3.46| )-( |3.48| ): 



{p,\x{t)x{t')\pa) = TjGUt.t'), 






{Pb I m Pit') I Pa ) = thG%{t, t') = ihG%{t\ t), 
{p,\p{t)p{t')\pa) =thMGUt,t'). 



(3.55) 



G%{t\t), (3.56) 



GUt\t). (3.57) 



(3.58) 
(3.59) 
(3.60) 



The relation between the similar-looking actions ( 3.43| ) and (|3.53|) becomes more transparent 
by reexpressing both in terms of partial derivatives of the classical solutions Xc\{t),Xc\{t), 
Pci{t) , Pc\{t) with respect to the end points Xb, Xa and Pb,Pa, respectively. In the configuration 
representation we obtain 



A{Xb,tb;Xa,ta)[k{t)J{t)] = -{Xb,Xa) 



I dpb_ dpb_ * 
dxb dxa 

dpa dpa 



\ dxb dxa J 



Xb 



yxaj 



tb 

+ dt (Xb, Xa) 



ta 



( dp,,{t) dx^ijt) \ 

dxb dxb 

dpdjt) dxdjt) 

V dXa dXa j 






1 dxh 



tb tb 



2dpa 



dt / dt'{k{t),j{t)) 



Ea ta 



e(t - f) 



dxa dxb dxa dxb 

dxci{t) dpciit') dxdit) dxdit') 



\ dxa dxb dxa dxb I 



+e(t' - 1) 



( dp,i{t)dp,iit') dpcijt) dx,i{t') \ 

dxb dxa dxb dxa 

dxd (t) dpc\ {t') dxc\ (t) dxc\ (f) 






(3.61) 



V dxb dxa dxb dxa j 

The momentum representation, on the other hand, has the analogues form with x and p 



interchanged: 



A{j)b,tb\Pa,ta){k(t),2(t)\ = - {pb,Pa) 



dxb dxb 

dpb dpa 

dxg dxg 

\ dpb dpa J 



^.^ 



\p.j 



•■b 
+ dt {pb,Pa) 



ta 



( dpait) dx^\ 

dpb dpb 

dpcijt) dxcijt) 

V dpa dpa J 



1 dxh 



tb tb 



2dpa 



dt / dt'{k{t),j{t)) 



ta ta 



( dp^iit) dp,i{t') dp,,{t)dx,,{t')\ 



e(t - 1') 



dpa dpb dpa dpb 

dxcijt) dpcijt') dxcijt) dxcijt') 
V dpa dpb dpa dpb j 



+e{t' - 1) 



( dpci{t)dpci{t') dpcijt) dxcijt') \ 

dpb dpa dpb dpa 

dxci (t) dpd (f) dxci (t) dxd it') 



^ Ht') ^ 



(3.62) 



V dpb dpa dpb dpa / 

These expressions for the generating functionals ( ^.42| ) and ( |3.50| ) exhibit clearly the sym- 
metry properties (Of) and (12] 



IV. SMEARING FORMULA FOR HARMONIC FLUCTUATIONS 



As a first application of the generating functional (|3.42|) we derive a general rule for 
calculating correlation functions of polynomial or nonpolynomial functions of x{t) and p{t). 
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The result will be expressed in the form of a smearing formula. This formula will represent an 
essential tool for calculating perturbation expansions with nonpolynomial interactions. Such 
expansions serve in variational perturbation theory to obtain convergent approximations for 
quantum-statistical partition functions |^ or density matrices 0. 

Consider the correlation functions of a product of local functions for vanishing currents 



( Fi(a;(ti)) F2ix{t2)) . . . F;v(a;(t,v)) FN+MtN+i)) F;v+2(p(t^+2)) • • • F^+A/(p(t7v+Af)) )n 



^b:'^o. 



Xl,,tl, 



Xl) 65 I Xa Zq^ j 



M 



VxVv ^ 

"^ n [Fn{x{tn))] n [FN+m{p{tM+m))] CXp ] ^^[p, X; 0, 0] L (4.1) 



^a i^-a 



27ih 



n=\ 



m=\ 



h 



where the harmonic time evolution amplitude with zero external currents {xf, tf, \xata)[0, 0] is 
written as (xf, tf, \ Xa ta)- By Fourier transforming the functions F„(x(t„)) and -Fiv+.m(p(iAr+m)) 

according to 

+00 +00 +00 

Fn{x{tn)) = J dXnFn{Xn)S{Xn-x{tn)) = J dx^ F{Xn) J — ^ CXp {i^„(x„ - x(r„))} 

—00 —00 —00 

(4.2) 



and 

FN+m{p(tN+m)) 

dp 



2Tih 



FN+miPm)5{Pm " PitN+m)) 



dPr, 
2TTh' 



FN+m{Pm) / c^/€^e-*'^'"(^™-^(*-+™«/^ (4.3) 



the correlation functions ( [4.1| ) may be reexpressed as 



{F,{x{h)) . . . FN+M{p{tN+M))yd 



^bi'^a 



1 ^ 

— - — n 

(, Xb i^b I -^a l^a ) n=l 



M 
rra=l 



dpn 
27Tn 



+00 



FN+m{Pm) / dK^e ^'^-P™/^ 



Ztt 



{Xbtb\Xata)[k,j], (4.4) 



where the generating functional is given by p.42|). The currents j{t) and k{t) are specialized 



to 
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N M 

jit) = -JlJ2 ^"^(^ ~ ^")' ^W = Yl '^rnSit - tN+m) 

n=l m=l 



(4.5) 



Inserting these equations into the action (|3.43 ) and the Green functions ( 3.22 ), ( 3.44 ) and 
( p.45| ), we find the Fourier decomposition of the the generating functional ( p.42|) , so that 
the correlation functions f|4.4|) become 



Fi{x{ti))... FN+MipitN+M)) )n 



'^hi-^a 



N 

n 

n=l 



+ 00 



/TT 



'^^" -i?n(a;„-a;d(tn)) 



M 

n 

m,=l 



X exp 



+00 



2Trh 



FN+miPm) I c?«:^e-*'^'"(^--^^'(*^+-))/^ 



■f. N N M -j.^ M 

/ . ^n^jj ^n' "•" ' Z^ Z^ ^n^jk ^m 1^ 2^ l^m^ kk 



2M Vi 

n,n'=l 



n=l m=l 



2h 



r\jYY)/ 



(4.6) 



rn,m'=l 



where we used the abbreviations 



G^j" — G^j{tn,tn'), G'^J!^ — G^j^{tn,tN^rn), G^fefc" — Glf.{tN^rn,tN+m') ■ (4.7) 



To proceed, it is more convenient to write expression ( |4.6| ) as a convolution integral 



^ r f+oa 1 *^ r r+°° dv 

F,{x{h))...FN+M{p{tN+M)))n""''=I[ / dXnFn{Xn) \[ / ^ FN+m{Pm 

nil L^-- J ^=1 K— 27r;i 

X ( -^ 1 P{xi,...,xn,Pi,...,Pm) (4.8) 



involving the Gaussian distribution 



P{xi,...,pm) 



(27r 



a V exp < «w V v G v 



lAT 



(4.9) 



The dimensionless vectors v and w have N + M components and are defined as 






IMQ 



IMQ 



MQ 



,N, 



1.1, ... , 



-I^M 



(4.10) 



and 
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w 



1 






VhMn 



xi - Xci(ti)] , • • • , \/ ^— [xn - Xdit 

- [Pl - Pcl{tN+l)] , • • • , 



N 



VfiMn 



{PM -Pc\{tN+M)) • 



(4.11) 



The {N + M)x {N + M)-matrix of Green functions 



< A B^ 



G 



\ 



B^ C 



(4.12) 



can be decomposed into block matrices A, B, and C. The A^ x A^-matrix A and the M x M- 
matrix C are defined by 



/ 



A = n 



jj jj 



(^12 ^11 
jj jj 



jj 



\ 



'^jj 



\ jj jj ' ' ' jj ] 



I 



^=a 



^kk ^kk 



^kk 



\ 



^kk ^kk 



ri2M 
^kk 



riiM r^2M r^u 

\ ^kk ^kk ' ' ' ^kk I 



(4.13) 



and yield quadratic forms of the position and momentum variables, respectively. The A^ x M- 



matrix 



B 



-'^jk 



r'21 
-'^jk 



r<i2 
-'^jk 



-'^jk 



y ^jk ^jk 



r<iM 
-'^jk 



r<2M 
''^jk 



'^jk j 



(4.14) 



gives rise to quadratic terms which are linear in both position and momentum variables. 



The multidimensional integral in ( [4. 91 ) is of Fresnel type and can easily be done, yielding 
and explicit expression for the Gaussian distribution 



P(Xi,... XAr,Pi,...,PM) 



1 



^iv+Af(27r)^-A^detG 



exp|-w^G"^w| , (4.15) 
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where G ^ represents the matrix inverse of ( |4.12|) whose block form is 

f ^ \ 



G-i 



with the abbreviation 



V 



X-1 -X-^BC-^ 

-C-^B^X-^ C7-1 + C-^B^X-^BC-^ 



(4.16) 



X = A- BC-'B' . 



llB^I 



Since the matrix G may be decomposed as 

X 
C-^B"" 1 
when the matrix C is regular, the determinant of G factorizes as follows 



G 



y,Cj 



\' 



J 



(4.17) 



(4.18) 



det G = det C det X. 
For singular matrix C but A regular, one may make use of another decomposition. 



(4.19) 



G 



V 



1 
S^A-1 X' 



\'ab^ 



(4.20) 



yO ly 



with X' = G — B A ^B. Then the determinant of G is given by 



det G = det X' det A. 



(4.21) 



With the Gaussian distribution (^4.15|) , our result ( |4.8|) constitutes a smearing formula v/hich 
describes the effect of harmonic fluctuations upon arbitrary products of functions of space 
and momentum variables at different times. 
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V. GENERALIZED WICK RULES AND FEYNMAN DIAGRAMS 

In applications, there often occur correlation functions for mixtures of nonpolynomial 
functions F{x{t)) or F{p{t)) and powers according to 

( F{x{t,)) x-{t,) )r% ( F{x{t,)) p-{t,) )r^ 
( FiPih)) x-{h) )-'-, ( F{p{t,))p-{h) )-•- . (5.1) 

In order to evaluate such correlation functions, we derive in this section generalized Wick 
rules and Feynman diagrams on the basis of the smearing formula (|4.8| ) . 

A. Ordinary Wick Rules 

It is well known that if one has to calculate expectation values of polynomials with even 
power, Wick's rule can be written as the sum over all possible permutations of products 
of two-point functions. We shortly recall to this expansion by considering the case of a 
position-dependent n-point correlation function, n even, defined as 

G(")(tl,...,t„) = (^(tl)■■■x(t„))?^-^^ (5.2) 

Note that it will be sufficient to study only the correlation functions involving the deviations 
from the classical path, respectively. This expectation value can be decomposed with the 
help of Wick's expansion 

pairs 

where p denotes the operation of pairwise index permutation. Thereby, the Green function 
^^^-'(^1,^2) is already given by ( |3.4(j| ). Note that Eq. ( [5. 3] ) may be considered as a consequence 
of a simple derivative rule 
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;F(x(ti))x(t2))r'' = {x{t^)i{t,)rr^{F'{x{t,))r^ 



-^fjj-^a 



(5.4) 



with F'{x) = dF{x)/dx. By applying this recursively, one eventually obtains ( |5.3| ). And 
conversely, the derivative rule ( p.4[ ) can be proved for j^o/ynomza/ functions F{x{t)), following 
directly from Wick's theorem (|5.3| ). 

The two-point Green functions G'*-^-'(ti, t2), occuring in (|5.3| ), can be considered as a Wick 
contraction which we introduce as follows: 



xiti)xih) = (5(ti)5(t2))c""" = '^G,,{hM). 

i{ti)p{t2) = {i{ti)p{t2)Y^'''^ =inGjk{ti,t2), 
I I 

P{ti)x{t2) = {p{ti)S:{t2))n"'''' =ifiGkj{ti,t2) = ihGjk{t2,h), 
I I 

P{ti)p{t2) = {Piti)pit2))n""'' = ihMGkkitiM)- 
I I 



(5.5) 
(5.6) 
(5.7) 
(5.8) 



Decomposing polynomial correlations of x{t) and p{t) with the help of these contractions 
corresponding to Eq. ( |5.3| ) or successively applying the derivative rule ( ^.4|) leads to following 
results 



;x"(t0:£™(t2))r" 

min(n,m) 

E Q 

Z=o,a+2,o+4,... 

:5"(ti)r(t2))r'' 

min(n,m) 

E Q 

Z=o,a+2,o+4,... 

;r(ti)5'"(t2))r" 



— Gj,,-(ti,tiJ 



{n-0/2 r 






— Gjj[t2,t2) 



-1 (m-0/2 



M 



Gjj{ti,ti) 



(n-l)/2 



\ihGjk{tiM)]' \i'^MGkk{t2M)] 



(m-l)/2 



min(n,m) 

^ Q [^/iMG'fc,(ti,ti)]('^-')/' [^/iG,fe(t2,ti)]' 

Z=a,a+2,o+4,... 



M 



Gjj{t2,t2) 



{m-l)/2 



'r{ti)p'^{t2)rr^ 



min(n,m) 

^ Q [^/iMGfcfc(ti,ti)]("-'^/' [^;iMGfcfc(ti,t2)]' [z/iMGfcfe(t2,t2)]^' 

=a,a+2,o+4,... 



(m-Z)/2 



(5.9) 



(5.10) 



(5.11) 



(5.12) 
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with the inultiphcity factor 



Q 



{n-l-l)\\{m-l-l)\\n\m\ 



(5.13) 



Note, that {—k)U = 1 for any positive integer k. For nonvanishing correlation, the sum 
n + m must be even so that the regulation parameter a is defined as follows: 



0, n,m even. 



a 



(5.14) 



1, n,m odd. 
The contractions defined in (|5.5|)-(|^) can be used to treat The desired derivative rules read 



F(5(ti))x"(t2))r 



E 



n\ 



l=a,a+2,a+4,... ^^ V--'- 

(F(5j(to)p"(t2))r" = 



— Gjj[t2,t2) 



(n~l)/2 



M 



Gjj{ti,t2) 



(F«(x(ti)))r", (5.15) 



E 



n\ 



l=a,a+2,a+4:,... V^ ij-.i- 

(F(p(to)r(t2))r'' = 



\ihMGkk{t2M)t-'^" \ihG,u{tiM)UF^'\5^{ti)) )T% (5.16) 



n\ 



E T^Tmn \inMG,,(t2M)t~'^" \ihMG,,{t,M)nF'^Hpiti)) )n"^", (5.17) 



n\ 



E 



— Gjj[t2,t2) 



-1 {n-l)/2 



l=a,a+2,a+4,... (^ '')■■''■ 

The parameter a distinguishes between even and odd power n: 



[^/^G,fc(t2,tl)]'(F«(p(tl)))^-^^ (5.18) 



a 



(5.19) 



0, n even, 

1, n odd, 

since even (odd) powers of n lead to even (odd) derivatives of the function F{x{ti)). The 
Ith derivative F^'-\x{ti)) is formed with respect to x(ti), and F^^\p{ti)) is the /th derivative 
with respect to p{ti). Note, that in the last line the Green function Gjk appears with 
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exchanged time arguments, which in this case happens to be inessential due to the symmetry 

Gjk{t2,tl) = Gkj{tl,t2). 

B. Generalized Wick Rule 



According to their derivation, the contractions ( 5.15 )-( 5.18| ) are only applicable to func- 
tions F{x[t)) and F{p{t)) which can be Taylor-expanded. In the following, we will show with 
the help of the smearing formula ( |4.8| ) that these derivative rules remain valid for functions 
F{x{t)) and F{p{t)) with Laurent expansions. Expectations of this type appear in varia- 
tional perturbation theory (see for position-position coupling Ref. [0]). Since the proceeding 



is similar in all the cases ( 5.15|) - (|5.18|) , we shall only discuss the expectation value 

(F(£(t0)p"(t2))r" (5.20) 

in detail. For this we consider the generating functional of all such expectation values 
following from ( [4.8|) 



(F(5(ti))e^P"(*^))^*'"" = ^^^ / dxF{x 



dp 
V- det G J -^-^^ '■^' J 2^ 



-e 



jp 



X exp 



Y^kkih, h) x^ - '2-Gjk{ti,t2) xp + -r-T--Gjj{ti,ti) p^ 



(5.21) 



.2detG' 
The p- integration can easily be done, leading to 

+ 00 

X / Fix + iTiG.Jh, h) j) e^*^"V2aG,,ft„t,) 

_i ^2'KihG,j{h,h)/M 

= e^^^^-(*-*^)^"/2 Y. T\ \^^G,,{UM)jUF^'\x{U)) )^''^^ (5.22) 
1=0 ''■ 

The correlation of two functions at different times has been reduced to a single-time expec- 
tation value of the Ith derivative of the function F{x(ti)) with respect to x(ti), denoted by 
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F^^\x(ti)), with Green functions describing the dependence on the second time. Expanding 
both sides in powers of j, we reobtain (|5.16|) . 

Now we demonstrate that the derivative rules ( ^.15|) - (|5.18|) for Laurent-expandable func- 



tions F{x{t)) and F{p{t)) also follow from generalized Wick rules. Without restriction of 
universality, we only consider the expectation value 

{F{S:{t^))x-{h)rr'^. (5.23) 

The proceeding to reduce the power of the polynomial at the expense of the function F{x(ti)) 

is as follows: 

la. If possible (n > 2), contract x{t2) x{t2) with multiplicity {n — 1), giving 

(n - 1) xih)xih) { F(£(ti)) x--\h) )^^'^% (5.24) 

else jump to lb. directly. 

lb. Contract F{x{ti)) x{t2) and let the remaining polynomial invariant. We define this 

contraction by the symbol 

F{x{t,)) ~x{t2) i-~\h) = x{h) x{t2) { F'{x{h)) x-\h) )r^ (5.25) 

I I I I I ^ ' 

Ic. Add the terms la. and lb. 

2. Repeat steps la.-lc. until only expectation values of F{x) or expectations of its deriva- 
tives remain. 

Summarizing, we can express the first power reduction by the generalized Wick rule 

{n>2) 



F{x{t,)) 5"(t2) )^^'^'' = {n-l) xih) xih) ( F(£(tO) x^^^h) f^^ 

I I 



+ F{5:{t,))x{t2)x-'\t2) (5.26) 
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with the contraction rules defined in (p.5|) and ( p.25| ). For n = 1, we obtain 

( F{x{t,)) x{h) )?^-- = mm ( F'{x{t,)) rr\ (5.27) 

which is vahd for any function F{x{t)) generahzing the rule (|5.4|) that was proved for poly- 



nomial functions only. Recursively applying this power reduction, we finally end up with 
the derivative rule ( ^.15 ). Note that the generalization of Wick's rule for mixed position- 



momentum or pure momentum couplings is done along similar lines, leading to the derivative 
rules (gl6D-(FJ^). 

C. New Feynman-like Rules for Nonpolynomial Interactions 

Higher-order perturbation expressions become usually complicated. For simple polyno- 
mial interactions, Feynman diagrams are a useful tool to classify perturbative contributions 
with the help of graphical rules. Here, we are going to set up analogous diagrammatic rules 
for perturbation expansions for nonpolynomial interactions V{x{t),p{t)), whose contribu- 
tions may be expressed as expectations values 

Jdtn---Jdt,{ Vix{Q,piQ) ■ ■ ■ Vixih), p{t,)) )^'""'" (5.28) 

From ( |5.5|) - (|5.8|) follows that we have four basic propagators whose graphical representation 
may be defined as (setting h = M = 1 from now on) 

h t2 = {x{h)x{t2))n""''' = iG^jitiM), 

t, t2 = {Pih) p{t2) )T^ = iGkkitiM), 

h — - t2 =(5(tl)p(t2))n""''=^G',fe(ti,t2), 

tl -— t2 = {Pih) X{t2)y^"'''' = lGkjiti,t2) = tGjkit2,ti). 
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A vertex is represented as usual by a small dot. The time variable is integrated over at a 
vertex in a perturbation expansion, 



dt. 

ta 



We now introduce the diagrammatic representations of the expectation value of arbitrary 
functions F{x{t)) or F{p{t)) and their derivatives as 

tb tb 

• ^ldt{F{i{t))r^-^^, ^ ^ jdt{F{p{t)))\ 

ta ta 

^b ^b 

i/ ^ \dt{F\i{t))r^^^^, i/ ^jdt{F'{p{t)))i 



\ ^bi-^a 

>n ! 



?'f^f + \\ \^b,^a 



ta ta 

^b ^b 

dt{F"{5:{t))rr\ <^ ^ jdt{F"{p{t))rr\ 

ta ta 



With these elements, we can compose Feynman graphs for two-point correlation functions 
of the type ( |5.1|) for arbitrary n by successively applying the generalized Wick rule (|5.26|) or 
directly using the derivative relations ( p.l5| )-( pTT8[ ). The general results become obvious by 



giving explicitly a graphical representation of the following four correlation functions 

th tb tb tb 

dh J dt2 ( F(£(ti)) £(t2) Vn""" = JdhJ dt2 2%(ti, ta) ( F\i{h)) )^'""'' (5.29) 

ta ^a ta ta 

= • • , 



tb tb tb tb ^ 

dt, I dh ( F{x{h)) x\t2) )^'""" = JdhJ dh I iG,,{t2, h) ( F(£(ti)) )T^ 

ta ta ta ta 

+ [.G,,(ti,t2)]^(F"(x(ti)))r"} (5.30) 

* <C:> ■ 

tb tb tb tb , 

rfti j dh { F(x(ti)) x'ih) rn""" = JdhJ dU 3 tG,,{h, h) tG,,ih, h) ( F'(5(ti)) )^^'"" 

ta ta ta ta 
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]3 / T?lll( ~(-L W \Xh,Xa 



+ [^G,,{t,,t,)r{F"\x{t,)))■r^^ (5.31) 



= 3 -A- 







dhj dh{F{x{h))x\t2)rn"'-= J dhj dhUtG,,it2,h)Y (F(5(ti)))?^^'^" 

+ [zG,,(ti,t2)]^(FW(5(tO))r4 (5.32) 




Mixed position-momentum and momentum-momentum-correlations and their graphical rep- 
resentations are given in Appendix |^. 

The consideration of higher-order correlations with more than one function F{x{t)) or 
F{p{t)) can be reduced to the results (|5.9| )- (|5.12| ) or (|5.15|) -( pT8D by expanding them with 
respect to the classical path or momentum, respectively. By expanding both functions in 
the expectation value, one obtains for example 



OO OO -I 

(Fi(5(t0)i^2(5(t2)))r" = E E ^7/W2,n(x'"(ti)x"(t2))r" (5.33) 



with 

h,m = F^^\0), z = l,2. (5.34) 

But, constructing graphical rules for such general correlations is more involved due to the 
various summations over products of powers of propagators Gjj{ti, tj) with i,j = 1, 2. 

Finally, we apply the diagrammatic rules to the anharmonic oscillator with x'^-interaction 
which is a powerful system being discussed in detail by the help of perturbation expansion 
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|1|, Chap. 3]. With the Green functions given by (p.22|) , ( p.44| ), and ( p.45|) , the two-point- 
correlation for anharmonic system with arbitrary time-dependent frequency can then be 
expressed graphically, yielding the known decomposition for the second-order perturbative 
contribution 



tb tb 



XfjjXa 



72 



24 




(5.35) 



with subscript c indicating that we restrict to connected graphs only. Beyond this, our theory 
allows to describe nonstandard systems with polynomial interactions ( ^.28|) depending on 
both, position and momentum, to higher-order. Finally, we want to give the graphs for a 
four-interation x^ p^ to second-order to see the variations of possible graphs in comparison 
with (|535|) : 

tb tb 

ta 



16 i 



> + 2 



+ 16 'I I I } +4 < 




. (5.36) 



We see, that we have the same class of graphs already occuring in (|5.35|) , however, with 
different propagators connecting the vertices. Thus, both classes decay into subclasses with 
different multiplicities, but the total numbers remain 72 and 24 for each type of class, 
respectively. Furthermore, all graphs are vacuum-like graphs. Eventually, it is easy to 
construct the Feynman graphs for polynomial correlations higher than second order by 
applying Wick's rule or the Feynman rules given in this section. 
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VI. SIMPLIFICATIONS FOR PERIODIC PATHS 

Up to now, we discussed the harmonic time evolution amphtude with arbitrary frequency 
and external sources j{t), k{t) and corresponding Green functions fulfilling Dirichlet bound- 
ary conditions. In the sense of the quantum mechanical partition function 



LLiAj \ Jb Lh tXj h rt 



(6.1) 



which is an integral over the time evolution amplitude for closed paths, it is of interest to 



investigate the generating functional for closed paths. In analogy to (|6.1| ), we define 



+ CO 

Z[j, k]= j dx{xtb\x to) [j, k] 



(6.2) 



with (|3.42|) for x^ = x^ = x. One immediatly observes that Z = Z[0,0]. The integral is 
easily done, giving 



Z[3M 



1 



th th 



^=^= exp < j dt I dt' 



I-a Cq, 



-mGut,t')m 



M-" 



+ j{t) G^^,{t,t') k{t') + k{t) Gl,{t,t')j{t') + Mk{t) Gt,{t,t') k{t') 



(6.3) 



The Green functions, expressed with fundamental solutions ( p.2|) , ( p.3|) , are found to be 



G^(t,t') 



1 
1 



G],{t,t') + 
Gkkitj't ) + 



a{ta,tb) 

1 

a{ta, tb) 

1 

a(ta,tb) 



9{t)g{t') 
9{t)m 

mgit'] 



Gkjit yt)y 



(6.4) 
(6.5) 
(6.6) 



with 



a{t,t') = Da{t')-Db{t)-2. 



(6.7) 
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Since the function 

git) = Da{t) + D,it) (6.8) 

is periodic, g(ta) = g(tb), due to conditions (|3.2|), (|373|), and ( p.6| ), also the Green function 
G^j{t,t') becomes periodic, 

G^^ita,t') = G^^iU,t'). (6.9) 

In analogy to the harmonic propagator without external sources ( [4.1D we can define expec- 
tation values consisting of A^ position-dependent functions and M momentum-dependent 
functions by 



( Fi(x(ti)) F2(x(t2)) ■ ■ ■ F^+MipitN+M)) h = 

1 r Dx 'Dv ( i ^ 

Zf ~M^ ^'^""^^'^^ F2{x{t2)) ■ ■■FM+M{p{tN+M)) exp |-^[p,x;0,0]| . (6.10) 



We remark that the generalization of Wick's rule and the graphical representation with the 
help of Feynman diagrams of such correlation functions is exactly the same as given in the 
last section after substituting the Green functions G{t, t') by G(t, t') and expectation values 

(O) by iwm . 



VII. SUMMARY AND OUTLOOK 

We have reduced generating functionals with fixed end points to those with vanishing 
end points by adding special singular sources to the currents. The new generating func- 
tionals were calculated explicitly for the harmonic oscillator with time-dependent frequency. 
From this expression, a smearing formula was derived which serves to calculate correlation 
functions for arbitrary polynomial or nonpolynomial position- and momentum-dependent 
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couplings. We have further found a generahzation of Wick's theorem of decomposing corre- 
lation functions involving functions of the canonic variables of the system. This gives rise 
to certain generalized Feynman rules for position- and momentum-dependent expectation 
values. 

Due to its universality, the theory should serve as a basis for investigating physical 
systems with nonstandard Hamiltonian via perturbation theory and its variational extension. 
Note, that a perturbation theory for momentum-dependent interactions arises in important 
field theories such as the nonlinear a-model. Our work is supposed to prepare the grounds 
for a more efficient perturbation treatment of such theory. 
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APPENDIX A: GENERALIZED CORRELATION FUNCTIONS 

In this appendix we give the expectations for the correlation between a general position 
or momentum dependent function and a polynomial up to order n = 4: 
Position-Momentum-Coupling: 

tb tb tb tb 

I dt, I dt2 ( F(£(ti)) Pit2) )T^ = JdtiJ dh tG,k{hM) ( F\x{t,)) )^-"» (Al) 



dt, I dt, ( F{S:{t^))p\t,) )^-"" = JdtiJ dt^ | tGkk{t2, h) ( F{x{t,)) )^-"'' 

la f-Q, ta ta 

l2 / I?"f^f+\\ \^b.Xa 



+ W,,{t,M)Y {F'\x{t,))Y^^-^\ (A2) 
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= * 



tb tb tf, tf, 




dh j dt2{F{5:{h))p\t2)Y^^''^ = j dh j dtA^iG,k{tiM)iGkk{t2M){F'{x{h))Y^'''^ 

+ {iGAti.h)f (F"'(Mti)))S-"] (A3) 

j dh j dh{F{x{h))p\hm^''- = j dh j dtA vGkk{hM)f (i^(£(ti)))n""" 

if i f ^ 

+ Q[tG,k{ti,h)f iGkk{t2M) (F"(£(ti)))n'""" 

+ VG,,{hM)f {F^'\x{t,))Y^A (A4) 



+ 6 W "i \ + ir'' ~> : 





Momentum-Position-Coupling: 

tb tb tb tb 

dh I dh { FiPih)) £(t2) fn'"'"' = JdtiJ dh tGk.ih, h) ( F\p{h)) )^-"" (A5) 

ta ia ta t-a 

= ^ , 

tb tb tb tb ^ 

dt, J dh { FiPit,)) 52(t2) )?^'""" = JdtiJ dh I ^G,,it2, h) { FiPit^)) )^-"" 

a ta ia ia 

+ [^G,Atuh)? {F"{p{U))yr^ (A6) 



tb tb tb tb ^ 

ta ia ia ia 

l3 / T^fff/-/, \\ \Xt-,.X 



+ VG,,{t,M)r {F"\p{t{))Tr^\ (A7) 
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= 3 






+ i^-—^- 









t-a t-a 






(AJ 



+ 6« 



«. 



Momentum-Momentum-Coupling: 



<6 <6 tb tb 



dt,jdt^{F{p{t,))p{t^)Y^ 



dh j dt2iGkk{tiM){F'(j>{h))Y^ 



(A9) 






dt,J dh{F{p{h))p'{h)r^ 



dh J dt2UGkk{t2,t2) {F{p{ti))y^ 

t t ^ 

+ [^G,,(ti,t2)]'(F"(p(ti)))r" 



Xi,,Xa 



(AlO) 




+ 




tb tb 

dt, I dt2{F{p{h))P\t2))\ 



th tb 



dt^ / dt2\3iGkkiti,t2)iGkkit2,t2){F'ipit^))) 



Xii^Xa 



ia ia 



3 / r^fff/-/ 



+ [tG,,ih,t2)r{F"\pit,))r^ 



\ Xfj^Xa 



(All) 



3 -i:^^^^ \ 




tb tb 

dt,jdt2{F{p{t,))p\t2)Y^ 

ta ta 



th th 



dh dt2\ [tGkkit2,t2)f (F(p(ti)))^""'^ 



ta ta 



6 [tG,k{h, t2)r ^Gkk{t2, t2) ( F"{p{U)) f^ 



+ WUhMt {F^'\P{h)))T^ 



^bi'^a 



(A12) 
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The case of position-position-coupling has aheady been calculated in Sect. |VC 
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